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ABSTRACT 

We present results of semi-analytic calculations which show clear evidence for changes in the non- 
equilibrium ionization behind a supernova remnant forward shock undergoing efficient diffusive shock 
acceleration (DSA). The efficient acceleration of particles (i.e., cosmic rays) lowers the shock temper- 
ature and raises the density of the shocked gas, thus altering the ionization state of the plasma in 
comparison to the test particle approximation where cosmic rays gain an insignificant fraction of the 
shock energy. The differences between the test particle and efficient acceleration cases are substan- 
tial and occur for both slow and fast temperature equilibration rates: in cases of higher acceleration 
efficiency, particular ion states are more populated at lower electron temperatures. We also present 
results which show that, in the efficient shock acceleration case, higher ionization fractions are reached 
noticeably closer to the shock front than in the test-particle case, clearly indicating that DSA may 
enhance thermal X-ray production. We attribute this to the higher postshock densities which lead to 
faster electron temperature equilibration and higher ionization rates. These spatial differences should 
be resolvable with current and future X-ray missions, and can be used as diagnostics in estimating 
the acceleration efficiency in cosmic-ray modified shocks. 

Subject headings: cosmic rays - thermal emission: ISM ~ shock waves - supernova remnants - X-rays: 
ISM 



1. INTRODUCTION 

In young supernova remnant (SNR) shocks, the ac- 
celeration of cosmic rays leads to a softening of the 
equation of state in the shocked plasma. This comes 
about because the diffusive shock acceleration (DSA) 
process turns some non-relativistic particles into rela- 
tivistic ones and because some of the highest energy rel- 
ativistic particles escape from the shock. Both of these 
effects lead to lower post-shock plasma temperatures as 
well as higher post-shock densit ies (e.g.. lJones fc EUisonl 
119911: iBerezhko fc Ellisonlll999[ ). The ionization state of 
shocked gas at a particular time is dependent upon both 
the gas density and the electron temperature. In light 
of this, DSA ought to leave its imprint on the ionization 
structure of the shocked gas. Toward this end, we present 
what we believe to be the first self-consistent model for 
SNR evolution which includes the hydrodynamics, the ef- 
fects of efficient shock acceleration, and a full treatment 
of the non-equilibrium ionization balance at the forward 
shock. 

A number of young SNRs show both nonthermal and 
thermal emission i n the r egion behind the forward shock, 
includi ng SN 1006 (iVink ct al. 2003; Baniba ct al. 200j, 
Tycho (iHwang et al.l 120021: JCassam-Chenai et aLlfeOOTn . 
and Kepler |Revnolds et al.l 2007t ). The thermal emission 



arises when the forward shock sweeps up the circumstel- 
lar medium (CSM) and heat s it to X-ray emittin g tem- 
peratures. As pointed out in lEUison et al.l (|2007ft (here- 
after DCE07), the thermal emission is often considerably 
fainter than the nonthermal emission, but there are cer- 
tainly examples where the thermal emission is as bright 
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or br ighter than any nonthermal emission (|Vink et al.l 
[2001) . In SNR RX J1713. 7-3946, the lack of thermal X- 
ray emission is an important constraint on the ambient 
density and significantly impacts m odels for TeV emis- 
sion (e.g., S lanc et al. 1999; Ellis on. Slane fc GaensleJ 
[200l iAliaronian~et aLll2007i : iKatz fc Waxmanll2008i r 

If the diffusive shock acceleration process in young 
SNRs is as efficient as generally believed, with > 50% 
of the shock ram kinetic energy going into relativistic 
particles, nonlinear DSA will influence the SNR hydrody- 
namics and be important for non-equilibrium ionization 
(NEI ) calculations (e.g., iDecourchelle. EUison fc Ballet) 
[2000HEnison fc Cassam-Chenaill2005D . DCE07 took the 
first steps in self-consistently coupling nonlinear DSA 
with NEI by tracking the electron temperature (Tg) and 
ionization age (defined as ngi, where rig is the electron 
density and t is the time since the material was shocked) 
as a function of time in hydrodynamic simulations of 
SNRs where the forward shock was efficiently produc- 
ing cosmic rays (CRs) and, as a result, was substan- 
tially modified from test-particle results. They found 
that, while both and net did differ between the test- 
particle and CR-modified cases, in the cases where DSA 
is highly efficient, the synchrotron emission in the X-ray 
range is considerably stronger than the thermal X-ray 
spectrum, and any differences in the thermal X-rays as 
a result of CR-modification are likely to be missed. In 
this paper, we extend the work of DCE07 by explicitly 
tracking the non-equilibrium ionization state in a CR- 
modified shock. The lower shock temperature and higher 
density that result from efficient DSA combine to shorten 
both the temperature equilibration and ionization equi- 
librium time-scale, and we show that this can have a 
dramatic effect on the ionization structure between the 
forward shock (FS) and the contact discontinuity (CD). 
Although we don't calculate the thermal X-ray emission 
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here, the cases we study show that efScient DSA can in- 
crease the ionization fraction of important elements and 
possibly enhance thermal X-ray emission. 

In § 2, we outline the changes to our model first pre- 
sented in DCE07 and discuss several caveats to our ap- 
proach. In § 3, we present our examples and discuss the 
quantitative and qualitative effects of efficient DSA on 
the ionization state and SNR structure. We also show 
how these effects might manifest themselves in current 
and future X-ray observations. In § 4, we summarize 
our results and outline our future enhancements to this 
model. 

2. CR-HYDRO + NEI MODEL 

Our spherically symmetric model uses the semi- 
analyt ic DS A calculation deve loped by lAmato fc Blasil 
(|2005f ) and iBlasi et al] (|2005( l and is similar to that 
used in DCE07, except that we now calculate the non- 
equilibrium ionization explicitly at every time step us- 
ing plasma parameters that are continually updated as 
the SNR evolves. In DCE07, the NEI was calculated at 
the end of the simulation using average plasma param- 
eters. We refer the reader to DCE07 for all details of 
the CR-hydro simulation apart for those discussed below 
detailing our dynamic NEI generalization. 

The DSA model us e d here differs from that describe d 
in lEUison et"an(|2007[ ). lEllison fc Cassam-Chenaj (|2005f ). 
and previous papers, in two important ways. First, we 
replace the "effective gamma," 7eff, approximation with 
a more realistic model of the effect escaping particles 
have on the shock dynamics. We now explicitly remove 
from the shocked plasma the energy that escaping par- 
ticles carry away from the forward shock. The ratio of 
specific heats of the shocked gas used in the simulations, 
7sk, is determined directly from the particle distribution 
function including the correct mix of rclativistic and non- 
relativistic particles. While the old effective gamma had 
the range 1 < -fcs < 5/3, the ratio of specific heats 7sk is 
constrained to lie between 4/3 and 5/3. These changes in 
the way escaping particles are treated, and 7eff is calcu- 
lated, become important for later stages of the SNR evo- 
lution, but do not produce significant cha nges in times 
as sho rt as 1000 yr. The results reported in lEUison et al.l 
()2007t l are not modified significantly by these changes. 

The second difference is that instead of specifiying a 
fi xed injection para meter, Xinj (this is ^ in equation (25) 
in lBlasi et al.|[2005[) . which then determines the acceler- 
ation efficiency, we now specify a fixed diffusive shock 
acceleration efficiency, edsa, and then determine Xinj ac- 
cordingly. This change makes the parameterization of 
the acceleration efficiency more transparent but does not 
change the basic approximation that is made. 

The semi-analytic DSA model we use does not cal- 
culate the acceleration efficiency self-consistently based 
upon the Mach number, the available acceleration 
time, and other relevant shock parameters; rather we 
parametrize the efficiency by Xinj, and the model then 
determines the shock structure self consistently. Fur- 
thermore, the DSA model assumes that the thermal par- 
ticles have a Maxwell-Boltzmann distribution with a su- 
pcrthcrmal tail. The actual shape of the quasi-thermal 
distribution, and the shape at the point where the su- 
perthermal tail joins it, are approximated since the semi- 
analytic calculation only self-consistently describes par- 



ticles with speeds greater than the shock speed, i.e., 
Vp ^ Wsk- The differences at low energies between what is 
assumed in the DSA model and the actual quasi-thermal 
distribution are expected to be small, but these differ- 
ences may become more important if the contributuion 
to ionization from superthermal particles is considered. 
Despite the approximations of the semi-analytic calcula- 
tion at quasi-thermal energies, it is the state-of-the-art 
since the actual quasi-thermal distribution can only be 
determined with plasma simulations and these are not 
yet available for SNR parameters. 

The ionization structure of shock heated gas at a par- 
ticular distance behind the shock in a SNR is determined 
by the electron density n,,, the electron temperature Te, 
and the ionization and recombination rates for each ion 
of interest. The structure is determined by solving the 
collisional ionization equations in a Lagrangian gas ele- 
ment behind the shock: 

lMp. = C{X^-\T,)fiX'-') + a{X\T,)f{X'+') 

Tie Dt 

-[CiX\n) + a{X'-\T,)]f{X^) . (1) 

Here, f{X^) is the fraction of element X in ion stage 
X"^ and G{X'\Tf,) and q;(A"', Tg) arc the ionization and 
recombination rates out of and into ion AT*, respectively. 

We calculate the electron temperature by assuming 
that the electrons ar e heated by C oulomb collisions with 
protons and helium (jSpitzeilll965[ ). We adopt this simple 
prescription, which gives a lower limit to the equilibra- 
tion time, knowing that the heating of electrons may, in 
fact, be far more complicated. For instance, there is rea- 
son to believe that coUisionless wave-particle interactions 
with the mag netic turbulence will be important (e.g., 
iLamingl l2001f ) , and recent work interpreting hydrogen 
line widths suggests that the electron-to-proton tempera- 
ture ratio behind some SNR blast waves depends mainly 
on the shock speed, a result implying a heating pro- 
cess substantially different from Coulomb coll isions (e.g., 
iGhavamian et al.ll2007t iRakowski et al.ll2008[ ). However, 
there remain large uncertainties in connecting the mea- 
sured line widths to the electron -to-proton temperature 
ratio (see iHeng fc SunvaevI l2008l ) . and until particle- in- 
cell (PIC) simulations are able to model non-relativistic, 
electron-proton shocks with parameters typical of SNRs, 
the plasr na physics of electron heating will r emain uncer- 
tain fsec lVladimirov. Bvkov fc El lison 20081 for a discus- 
sion of the limitations of PIC simulations in this regard). 
In order to model some of the complexity of electron 
heating, we scale the Coulomb equilibration time with a 
parameter, /cq, defined in Eq. ([3]) below. 

At the start of the simulation, we assume that the un- 
shocked electrons and ions are in equilibrium at a tem- 
perature To = 10^ K. We also assume that unshocked 
H and He are both 10% singly ionized and all heavier 
elements are initially neutral. While we note that this 
is not the precise equilibrium ionization state for 10^ K, 
we emphasize that none of our results depend in any 
significant way on the ionization state of the unshocked 
material as long as it is not fully neutral. In all of the 
results shown here we fix the helium number density at 
10% of the proton number density, Up^. 

At each time-step, we track the ionic state A' within 
each spherically symmetric fluid element by solving the 
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timc-dcpcndcnt ionization equations for each abundant 
element (H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe, and 
Ni). We solve the co upled set of equations wit h atomic 
data extract ed from iRaymond fc SmithI ()1977[ ). as first 
presen ted in lGaetz et al.l ( 1988f ) and updated by lEdgarl 

In Figure [T] we show an example of the time evolu- 
tion of the ionization fraction, /(X*), of high ionization 
states of oxygen (0^+, 0^+ and 0®+) in a mass shell that 
is crossed by the forward shock 100 yr after the explo- 
sion. For this example, as in all we show in this paper, 
we have fixed parameters typical of Type la supernovae, 
i.e., the kinetic energy in ejecta from the supernova ex- 
plosion EsN = 10^^ erg, the mass of the ejecta Mej = 
1.4M0, the density of the ejecta follows an exponential 
densit y profile as is gene rally assumed for Type la super- 
novae (|Dwarkadasll200G[ ). and we assume the supernova 
explodes in a circumstellar medium (CSM) which is uni- 
form with proton number density Jip.o and magnetic field 
strength Bq.^ In all of the models shown here, we take 
Bq = 15 /xG.^ The figure shows that the density and pro- 
ton temperature in the shell arc dropping with time as 
the electron temperature increases due to Coulomb col- 
lisions. After 1000 yr, the material is close to ionization 
equilibrium for these ions. 

Figure [T] also compares results for test-particle (TP) 
and efficient DSA. In all of the examples in this paper, 
we define TP acceleration as being 1% efficient, i.e., 1% 
of the ram kinetic energy of the forward shock is placed 
into superthermal particles. For all of our efficient ac- 
celeration cases, we assume 75% of the shock ram ki- 
netic energy is placed into superthermal particles, i.e., 
cdsa = 75%. Figure [1] shows that efficient DSA produces 
a higher postshock density and lower postshock temper- 
ature, as expected. What is also clear is that the high 
ionization states of oxygen become populated sooner in 
the eDSA = 75% case. This implies that, instead of sup- 
press ing thermal X-ray emission as has been suggeste d 
(e.g., iDrurv et al.ll2008t iMorlino. Amato fc Blasil 120081 ). 
efficient DSA can possibly enhance it. 

We make the following approximations in the NEI cal- 
culation, noting that these are in addition to approx- 
imations m a,de in the underly ing CR-hydro model (as 
described in lEUison et al.|[2007l . and references therein): 



• We assume that only electrons from the thermal 
population contribute to the non-equilibrium ion- 
ization. In nonlinear DSA, the energetic popu- 
lation emerges smoothly from the thermal pop- 
ulation (a nice examp le from a relativ istic PIC 
simulation is given in iSpitkovskvl [200l and su- 
pert hermal particles ma y contribute to ionization 
(see iPorguet et al. ^ '2001', for a tcst-particlc calcu- 
lation involving a MaxwcU-Boltzmann distribution 
with nonthermal tail). As we discussed above, su- 
perthermal particles are expected to contribute to 



^ We refer the reader to lEllison et al.l II2007I ) for a full discussion 
of the additional parameters required for the CR-hydro model. 

* This value for Bq is somewhat higher than the typically as- 
sumed 3 fiG and reflects the possibility that magnetic field ampli- 
fication (MFA) may be taking place. We emphasize, however, that 
we do not include MFA in the DSA calculation performed here. A 
large upstream magn etic field, Bp, will reduce th e effects of efficient 
DSA, as described in lBerezhko fc EUisonl llT999h . 



the ionization at some level. However, the sig- 
nificance of this nonthermal ionization, in shocks 
undergoing efficient particle acceleration, has not 
yet been determined and remains an area of active 
work. For the purposes of this paper, we assume 
any nonthermal contribution is small. 

• We only model the interaction region between 
the forward shock and the contact discontinuity 
where we assume cosmic elemental abundances. 
One reason for emphasizing the forward shock is 
that it is not certain that significant CR produc- 
tion occurs at the reverse shock i n SNRs (e.g., 
lEUison. Decourchelle fc Balletll2005[ ). 

• We only consider young SNRs and do not include 
the effects of radiative cooling. In the high-density 
limit, radiative losses could be significant and the 
cooling time-scale could be comparable to other dy- 
namical time-scales. We will investigate this effect 
in a subsequent paper. 

3. RESULTS 

In the following examples we investigate the effect the 
acceleration efficiency, edsa, and the CSM proton den- 
sity, Upfi^ has on the non-equilibrium ionization state of 
some selected elements. 

3.1. Ionization vs. Position 

In Figure [21 we plot the ionization fractions of 0^+ 
and 0^+ and Si^^+ and Si^^+ in the top two panels as a 
function of position behind the forward shock (FS). In all 
panels, test-particle results (eosA ~ 1%) are shown with 
dashed curves and efficient DSA results (edsa = 75%) 
are shown with solid curves. The electron density and 
electron and ion temperatures are shown in the bottom 
two panels.^ 

As the top two panels clearly show, higher ionization 
fractions are attained closer to the shock front in the ef- 
ficient DSA cases, as compared to the TP cases. For in- 
stance, in the efficient case, the fraction of 0^+ peaks at a 
distance R/Rps — 0.98 behind the shock, while in the TP 
case, this fraction peaks at R/Rps — 0.97. We attribute 
the increased ionization fractions closer to the shock as a 
direct result of higher postshock densities in the efficient 
DSA case. Note that the curves extend from the forward 
shock back to the contact discontinuity, indicating that 
the region between the forward shock and contact discon- 
tinuity is considerably narrower in the efficient accelera- 
tion case. This effe ct produces important morphological 



Decourchelle. Ellison fc Balletl 120001: 



Cassam-Chenai et al.ll2008f ). 



consequences (e.g.. 
IWarren et ani2005i : 

In Figure [31 we show the same quanitities as in Fig- 
ure [2l except that rip^o = 0-1 cm~'^. The lower CSM 
density results in lower shocked densities and in less 
rapid coUisional ionization behind the FS. For the ions 
we show, higher ionization states (i.e., 0*"+ and Si^^+) 
are considerably less populated downstream from the FS 
when Up^o is small. The differences resulting from DSA 
are less prominent but still evident; e.g., with = 0.1 
cm""^, 0^+ peaks behind the shock at R/Rps ~ 0.98 for 
the efficient case, and at ~ 0.95 in the test particle case. 

^ In all results shown, we assume that shocked protons and other 
ions have the same temperature. 
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To emphasize the importance of the different spatial 
structures of ionization with eosA and rip^o, we show, in 
Figure IH a closeup view of the shock fronts in Figures [5] 
and [3l Here, we have plotted the ionization fractions 
as functions of angular distance behind the shock, as- 
suming a distance of 1 kpc. In the high density case 
{np.o — Icm"'^; top panel), the fraction of O^"*" peaks 
right behind the shock at ~ 2" downstream, while it 
peaks ~ 5" behind the shock in the test particle case. 
In the lower density case (?T.p,o = 0.1 cm~'^; lower panel), 
0^+ peaks ~ 30" behind the shock in the efficient case, 
but peaks well beyond 50" behind the shock in the test 
particle case. Similar results are found for silicon. While 
these models are not scaled to match any particular 
Galactic SNR, we believe the angular separations shown 
here would be easily resolvable in current and future 
space-based X-ray observatories even when line-of-sight 
effects are taken into account. Thus, measuring the rela- 
tive fraction of H-like, He-like, and even Li-like charge 
states would provide a useful diagnostic in studies of 
Galactic SNRs undergoing efficient shock acceleration. 

Another interesting feature seen in Figures [2] and [3l 
is that the electron temperature is almost independent 
of eosA and only varies by a factor of ~ 2 between the 
n-pfi = 1 cm~'^ and Hp^Q = 0.1 cm~^ cases. This is in 
contrast to the ion temperatures, where generally lower 
ion temperatures occur in the higher density models, 
due to the lower shock Mach number, and where the 
large eosA cases have considerably lower ion tempera- 
tures than the test-particle cases. The fact that lower 
postshock te mperatures o ccur in efficient DSA is well 
known (e.g.. lEllisonll2000l ). The electron temperature is 
influenced by this and by the higher densities that occur 
with efficient DSA. The higher postshock densities im- 
ply more collisions between electrons and ions, and thus 
more rapid temperature equilibration. The higher elec- 
tron temperature combined with the higher postshock 
density leads to more rapid ionization, and thus higher 
charge states closer to the forward shock. 

3.2. Ionization vs Equilibration Timescale 

As is clear from Figures[Tl[2l and[3l the ionization frac- 
tion for high charge state ions can increase with acceler- 
ation efficiency. Since the electron temperature is almost 
independent of rip o in these cases, we attribute this effect 
mainly to the higher postshock densities. However, we 
have assumed a particular model for temperature equili- 
bration between protons and electrons, namely that elec- 
trons start off cold and equilbration with the hot protons 
occurs only through Coulomb collisions whe re the equi- 
libration timescale is given by (lSpitzej [l965l Eq. 5-31): 

, _ irnpUiekJ / Tp^ \ ^ (o) 

~ 8(27r)i/2npZ2Z2e4lnA \mp me J ' ^ ' 
Here, rup is the proton mass and Tp is the shocked pro- 
ton temperature and d efinitions of the other terms are 
given in ISpitzeil (|1965f ). It's important to note that 
Eq. ([2]) places strict limits on how low the electron to 
proton temperature ratio can be behind the shock (see 
I Hughes. Rakowski fc Decourchellell2000() ; if other equili- 
bration mechanisms are important, such as plasma wave 
interactions, equilibration will occur more rapidly. To 
investigate the effects of more rapid temperature equili- 
bration, we define a parameter. < fcq < 1, and use the 



equilibration time t'^^ in our calculations where, 

^eq ~ fcqtcq ■ (3) 

In the results shown in Figures [1] [21 El and [H we have 
assumed fcq = 1. 

In Figure O we compare the ionization fraction of 
0^+ for edsa — 1% and edsa = 75% calculated with 
fcq ~ 1 (black curves in all panels) and f^q ==0.1 (red 
curves in all panels). For both values of edsa, f{0'^~^) 
is larger immediately behind the shock for rapid equili- 
bration (fcq = 0.1) but drops below the /eq = 1 value 
further downstream as 0®+ becomes populated. The 
temperature plots in the bottom two panels show that 
the electrons and protons have come into equilibrium 
for a range of radii (i.e., 0.86 < R/Rps S 0.98) when 
^dsa — 75% and fcq ==0.1, but remain far from equilib- 
rium for fcq = 1 regardless of edsa- The equilibration 
rate changes the ionization structure for this particular 
ion, producing changes comparible in scale to those pro- 
duced by efficient DSA. 

To quantify these effects further, we look at a point 
midway between the contact discontinuity and FS, i.e., at 
R/Rfs ^ 0.89 for cusa = 75% and at R/Rfs ^ 0.83 for 
cdsa = 1% in Figure [HI At these locations, the electron 
to proton temperature ratios are: {Te/Tp)'Fp ~ 0.11 and 
iTe/Tp)^L ^ 0.36, for fcq = 1, and {TjTp)^^ ^ 0.3 
and (T'e/rp)NL = 1 for fcq = 0.1, i.e., the ratios are 
about 3 times larger with rapid equilibration. At these 
midpoint locations, the ionization fractions of O''^ range 
from /(0^+)~ 0.05 for fcq = 1 and edsa = 75%, to 
/(0^+)~ 0.23 for fcq = 0.1 and edsa = 1%, i.e., about 
a factor of five span. 

The electron temperature ratio for fcq = 1 is 
(re,NL/re,Tp)/e,=i = 1 .8 X lO" A72.5 X lO" X ^ 0.7 
and the ratio for fcq = 0.1 is (Te^NL/T'e,Tp)/eq=o.i = 
ixlQ^ K/QxlQ"^ K ~ 0.5. For the particular parame- 
ters used in this example, the electron temperature stays 
within a factor of ~ 2 for a wide spread in edsa and equi- 
libration time, while /(0''+) varies by a factor of ~ 5. 

3.3. Emission Measure vs. Acceleration Efficiency 

As seen in Figures[2]or[31 the plasma density is greatest 
immediately behind the shock where the electron tem- 
perature is lowest. Since the rate for electron temper- 
ature equilibration depends on the proton temperature 
and density and both the temperature and density de- 
pend on CDS a, the NEI calculation will depend in a com- 
plicated fashion on the forward shock dynamics and the 
evolution of the interaction region between the CD and 
FS. Of course, the important property is the emission 
the plasma produces and this can be characterized by 
the emission measure (EM) and the differential emission 
measure (DEM). 

In Figure [51 we plot the emission measure for individ- 
ual ions, EM = NxfiX\R)nenpdV, and in Figure [3 
we plot ionic differential emission measures, DEM — 
J2 Nxf{X\ R)nenpdV/d{\ogTe), where Nx is the abun- 
dance of element X relative to hydrogen, /(AT', R) is the 
ionization fraction for the ion X* at a distance R behind 
the shock, and dV is the volume of the shell where EM 
or DEM is determined. The EM plotted in Figure [SI is a 
line-of-sight projection normalized to 1 cm^ surface area. 
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and the DEM is obtained by summing over the region 
between the CD and FS. 

Figure [S] clearly shows that the emission for these ions 
peaks much closer to the FS and is considerably stronger 
with efficient DSA than in the TP case. Figure [7] shows 
that the peak emission for these two ions shifts down in 
temperature by about a factor of 2 ('^ 1 keV) when 
efficient DSA occurs. These two effects are quite signifi- 
cant for individual ions and should be observable. Nev- 
ertheless, the emission from a full set of ions needs to 
be calculated and the results folded through a detectors' 
response before the signature of efficient DSA can be 
quantitatively determined. 

4. DISCUSSION AND CONCLUSIONS 

We have presented a calculation of non-equilibrium 
ionization in a hydrodynamic simulation of SNRs un- 
dergoing efficient DSA. While we have only explored a 
limited range of parameters in this paper, it's clear that 
the production of CRs by the outer blast wave modifies 
the SNR evolution and structure enough to produce sig- 
nificant changes in the ionization of the shocked material 
between the forward shock and contact discontinuity. In 
particular, higher ionization states are reached at lower 
electron temperatures (compared to the test particle case) 
because of the increase in post shock density due to the 
increased shock compression. The calculation of thermal 
X-ray line emission requires the additional step of cou- 
pling the resultant ionization state vectors to a plasma 
emissivity code, work which is in progress. Nevertheless, 
our results clearly show that taking DSA into account 
and dynamically calculating the NEI produces changes 
in the ionization fractions of important elements that 
should translate into noticeable changes in the interpre- 
tation of X-ray line emission observed from young SNRs. 

Our main results are the following: 

• Compared to the test-particle case, the increase 
in ionization that accompanied DSA in our ex- 
amples suggests that efficient DSA will result in 
an increase in the overall thermal X-ray emission 
(see Figure [H]) . We note that an increase in ther- 
mal emission with increasing acceleration efficiency 
is evident in our earlier results which explored a 
slightly different p arameter space (i.e., Figures 7 
and 8, lElhson et a l. 2007,). The actual increase 
may depend importantly on other model param- 
eters, such as the CSM density, and it is important 
to explore a more expanded parameter space to de- 
termine how broadly valid our results are. This 
work is in progress. However, regardless of whether 
or not efficient DSA increases the integrated ther- 
mal emission over the test-particle case, some ther- 
mal emission is expected because ionization is not 
suppressed when efficient DSA occurs. As Fig- 
ure [T] shows, electrons reach X-ray emitting tem- 
peratures well before they come into equilibration 
with protons and nearly as rapidly with or without 
efficient DSA. This occurs even if only Coulomb 
equlibration is ass umed. This is in contrast to 
recent clainis (e.g ., iMorlino. Amato fc Blasil 120081 : 
iDrurv et al.l 12008!) that very weak thermal X-ray 
emission might result from efficient shock acceler- 
ation. 



• Compared to the test-particle case, ionization oc- 
curs more rapidly and, therefore, closer to the FS, 
with efficient acceleration (see Figures 2] and 
The differences in spatial structure should be large 
enough to observe and may be used as a discrimi- 
nant for the level of CR-modification, if a particular 
ion state is coupled to other known properties, such 
as the dynamics and ambient conditions. 

• Efficient DSA leads to more efficient Coulomb heat- 
ing of electrons and faster equilibration with ions, 
relative to the test particle case. This results be- 
cause the shocked plasma temperature is lower and 
the shocked density is higher when efficient DSA 
occurs. We showed, with a simple parameteriza- 
tion of the thermal equilibration time, that the sig- 
nature of efficient DSA on the ionization state re- 
mains apparent for equilibration more rapid than 
occurs with just Coulomb collisions. 

• Using the differential emission measure, we showed 
that the maximum emission from a particular ion 
state occurs at a significantly lower electron tem- 
perature with efficient DSA. For the ions shown in 
Figure [71 the difference in Tg for peak emission is 
on the order of 1 keV while the maximum DEM re- 
mains almost constant. A difference this large will 
have an important impact on the interpretation of 
thermal X-ray emission in young SNRs. 

Currently, we do not treat radiative or slow shocks, 
but these regimes are easily explored. For instance in a 
radiative shock, the cooling time might be comparable 
to the energy loss time in a cosmic ray modified shock. 
Increases in the density will enhance the cooling to the 
point where radiative losses mig ht rival losses from effi- 
cient DSA (jWagner et all 120061 ). We intend to explore 
this regime in a forthcoming paper. 

While we only considered shocked CSM here, we will 
consider shocked ejecta in future work. In the ejecta, 
the electron density can be higher and the tempera- 
ture may be lower but, more importantly, the abun- 
dance structure is far more complicated than for CSM 
and calculations of X-ray emission are intrinsically more 
difficult. Furthermore, simple arguments based on the 
expansion of the ejecta material suggest that the mag- 
netic field may be too low to support DSA by the reverse 
shock. Nevertheless, there has b een speculation that 
particles are accelera ted there (e.g., [G otthclf et al.J 20011 : 
lUchivama fc Aharonian, 20081 iHeldeT fc V ink 200a ) and 
if DSA is efficient at the reverse shock, it will likely alter 
the ionization balance of the shocked ejecta as much as 
shown here for the shocked CSM. 

Finally, while we have limited our examples here to 
SNRs expanding into a uniform medium typical of Type 
la supernovae, we emphasize that a wider parameter 
space should be explored, in terms of both the struc- 
ture of the ambient medium (i.e. pre-SN winds) and the 
parameters which determine the cosmic ray acceleration 
efficiency. These cases will be addressed in a follow-up 
paper. 
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Age [yr] 

Fig. 1. — Time evolution of a spherically symmetric Lagrangian mass shell which is crossed by the forward shock at 100 yr. The top 
panel shows the evolution of high ionization states of oxygen, the middle panel shows the electron number density, and the bottom panel 
shows the electron and proton temperatures, assuming Coulomb equilibration. In all panels, the solid curves correspond to a model with 
75% DSA efficiency, while the dashed curves are for a TP model with eosA = 1%. The CSM proton number density for this example 
is rip = lcm~^. Here, and in all other examples, the unshockod CSM temperature is Tq = 10 K, and the unshocked magnetic field is 
Bo = 15^lG. 
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Fig. 2. — Spatial profiles of H- and He-like oxygen and silicon, electron density, and temperature as a function of distance behind the 
forward shock. In the bottom panel, the curves labeled Ti are ion (or proton) temperatures and those Te arc electron temperatures. Here, 
and in figures l3l5l that follow, we show values from spherically symmetric shells as a function of R or AR, not line-of-sight projections. In 
all panels, solid curves correspond to models with 75% efficiency, while the dashed lines correspond to TP models. These models arc for a 
CSM proton density of rip^o = 1 cm~^ and are calculated at tsNR = 1000 yr. In the model with 75% efficiency, the forward shock velocity 
is ~ 1800 km s~^, while in the test particle model, it is ~ 2200 km s~^ at tsNR = 1000 yr. 
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Fig. 3. — Spatial profiles of oxygen and silicon ions, electron density, and temperature as a function of distance behind the forward shock. 
In all panels, solid curves correspond to models with 75% cfhcioncy, while the dashed lines correspond to TP models. These models are for 
a CSM proton density Upfi = 0.1 cm"'' and are calculated at tsNR = 1000 yr. In the model with 75% efficiency, the forward shock velocity 
is 3200 km s~^, while in the test particle model, it is ~ 3600 km s~^ at tgNR = 1000 yr. 
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Fig. 4. — Top: Ionization fraction as a function of distance behind the forward shock for O®"^ and O^"*" with np,o = 1.0cm ^ . Bottom: 
Ionization fractions of O^"^ and O^^" with Jip^o = 0.1 cm~'^. In both panels, the solid curves are for edsa = 75% and the dashed curves are 
for edsa = 1%- The angular scale is determined assuming the SNR is at a distance of 1 kpc and the results are calculated at igt-ja = 1000 yr. 



Non-Equilibrium Ionization in Cosmic Ray Modified Shocks 



11 



0.8 



X 



-*-> 

u 
k_ 

CL 

E 

0) 



0.6 - 

0.4 - 

0.2 - 


0.6 - 

0.4 - 

0.2 - 



10° 



1 — 1 1 n7 



10' =- 



lO** 



10° r 



10' r 



T" 



T" 



«DSA="75% black: f,„=1. red: f.„=0.1 



T" 



eq 



^7+ 




I I I I I I I I I I I I 

black: f =1. red: f =0.1 



eq 




I I I I I I I I I I I I 



I I I I 



solid: T , dashed: T 




Fig. 5. — Ionization fraction and temperature calculated between the contact discontinuity and FS. All calculations are at tsNR = 1000 yr 
and assume rip o = lcm~^. In all panels, black curves assume /eq = 1 and red curves assume /cq = 0.1. In the bottom two panels, the 
solid curves are the shocked electron temperature, Te, and the dashed cuvres arc the shocked proton temperature, Tp. As in Figures [2l and 
131 the left end of each curve is at the position of the CD. 
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Fig. 6. — Line-of-sight projection of the emission measure (EM) for O^"^ and Si^'^^ as labeled. The solid curves are for edsa = 75% and 
the dashed curves are for edsa = !%• The angular distance, AR, from the FS is determined assuming the SNR is at 1 kpc and the results 
are calculated at tsNR = 1000 yr with rip.o = Icm""^ and icq = 1. 
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Fig. 7. — Differential emission measure (DEM) vs. electron temperature for and Si^^"*" as labeled. The solid curves are for 

£dsa = 75% and the dashed curves are for eoSA = 1%. The results are calculated at tsNR = 1000 yr with ripfi = 1 cm~^ and teq = 1. 



